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Output-error parameter estimation, normally a post-flight batch technique, was applied 
to real-time dynamic modeling problems. Variations on the traditional algorithm were 
investigated with the goal of making the method suitable for operation in real time. Im- 
plementation recommendations are given that are dependent on the modeling problem of 
interest. Application to flight test data showed that accurate parameter estimates and un- 
certainties for the short-period dynamics model were available every 2 s using time domain 
data, or every 3 s using frequency domain data. The data compatibility problem was also 
solved in real time, providing corrected sensor measurements every 4 s. If uncertainty 
corrections for colored residuals are omitted, this rate can be increased to every 0.5 s. 


Nomenclature 


Roman 


A 

perturbation value 

A,B,C,D 

linear system matrices 

A 

elevator deflection [rad] 

1 

body translational accelerations [ft/s^] 

a 

standard deviation 

cov 

covariance 

4>, 6, V” 

Euler angles [rad] 

9 

gravitational acceleration [ft/s^] 

CO 

frequency [rad/s] 

I 

identity matrix 



j 

imaginary number, \/—l 

Superscripts 

P,9N 

body roll, pitch, and yaw rates [rad/s] 

-1 

matrix inverse 

q 

dynamic pressure [Ibf/ft^] 

T 

transpose 

t 

time [s] 


time derivative 

U, V, w 

body translational velocities [ft/s] 

t 

complex conjugate transpose 

V 

airspeed [ft/s] 


estimated value 

Greek 


Subscripts 


a 

angle of attack [rad] 

0 

trim or starting value 

p 

sideslip angle [rad] 



7 

flight path angle [rad] 
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I. Introduction 


O UTPUT error has matured into one of the most widely accepted and used methods in the aerospace 
industry for estimating unknown model parameters from measured flight data. It is a straightforward 
method based on maximum likelihood estimation that uses realistic assumptions consistent with standard 
flight test practices and the real physical situation, and also possesses useful analytic properties. These fa- 
vorable characteristics have resulted in successful application to a broad range of research problems including 
aerodynamic modeling, control law verification and validation, envelope expansion, flight path reconstruc- 
tion, and accident investigation. Popularization of this method is due, in part, to numerous researchers 
releasing software packages that combine the rigorous estimation theory with pragmatic numerical methods 
and practical engineering experience. 

Aircraft system identification is currently undergoing a paradigm shift towards real-time modeling during 
flight. This effort is facilitated by improvements in hardware and computing capability, and driven by the 
need to reduce costs and improve safety. For example, immediate feedback of modeling quality during control 
law flight tests can lead to shorter maneuver lengths, or can diagnose problems while on condition, thereby 
saving time and costs while avoiding subsequent tests. Safety issues, such as icing or fault-detection, can be 
addressed in real time by monitoring estimation results.® 

There are however several aspects of output error that make a real-time implementation difficult. Output 
error is a parametric modeling technique that requires prior knowledge of a model structure. This requirement 
is problematic if there is little or no information before a flight test or if the dynamics change, for example 
because of structural damage to the aircraft. Another difficulty is that output error is a nonlinear estimation 
problem that uses an iterative multi-dimensional nonlinear optimization method to determine the solution. 
These iterative methods not only require more computation time than non-iterative solutions, but they 
require good starting values for the unknown model parameters. Poor starting values can lead to excessive 
iterations, incorrect solutions, or divergence. Another problem is that after a solution is found, additional 
effort is needed to correct uncertainty estimates for colored residuals, which appear in practice because of 
unmodeled dynamics but are not assumed in the theory.® A. 8 penally, it is exceedingly difficult to apply 
traditional output error to unstable aircraft, such as fighters or damaged aircraft, because small errors in the 
model parameters and initial conditions quickly lead to divergent time histories of the model outputs.®’® 
The closest attempt to a real-time implementation of output error was by Hui et al.,^® who used FOR- 
TRAN software onboard a Dassault Falcon 20 to produce parameter estimates in “near real time,” soon 
after a maneuver completed using a batch process. Other researchers have looked to alternative methods for 
real-time parameter estimation. Filter-error and the extended Kalman filter have been used, especially when 
turbulence is present, but involve tuning parameters and produce inaccurate uncertainties. Equation- 
error parameter estimation in the frequency domain has been applied successfully in real time,^® but this 
method requires an extra step in the time domain to estimate bias parameters, and does not directly opti- 
mize the model parameters to fit measured output data in the time domain. A real-time frequency response 
method can be used without knowing the model structure, but additional effort is needed to deduce the 
structure from the results and estimate parameters. Currently there is no working real-time parameter 
estimation technique that uses the output-error approach. 

This paper addresses some of the shortcomings of output error, to make it suitable for solving real-time 
parameter estimation problems. Variations to the method are implemented and tested to determine best 
practices, and then results are applied to flight test data to estimate model parameters in real time. Section II 
presents traditional output error and its variations. A subscale flight test vehicle and its nonlinear flight 
simulation are described in Section III, as well as a maneuver used for system identification. Output-error 
results are presented in Section IV for the real-time estimation of short period model parameters and for 
solving the data compatibility problem during flight. Discussion is provided in Section V. 

II. Methods 


A. Traditional Output Error 

This section introduces traditional output error. For a comprehensive description, see Refs. [5, 15]. 
Consider the continuous-time model of the aircraft flight dynamics 

x(t) = f (x(t),u(f),0) , x(fo)=xo (1) 


2 of 16 


American Institute of Aeronautics and Astronautics 



where x is the state vector, u is the applied control input, and 0 is a vector of Up unknown model parameters. 
Disturbances in the state equations, e.g. due to turbulence or model structure errors, are neglected, and 
the control input is assumed to be measured without error. Consider also the no discrete-time sensor 
measurements 

z(*) = y(i) + i = l,2,...,N (2) 

observed at N regular time intervals, where y is the model output and is a Gaussian white noise sequence 
with zero mean and covariance R. 

The objective of output error is to determine the unknown model parameters that best match the model 
output to the measured output. This is done by minimizing the cost function 

1 ^ 

where 

■u(t) = z{i) - y(i) (4) 

is the modeling residual. Optimization requires an iterative method because the relationship between the 
model parameters and the model outputs is nonlinear. This process is outlined as follows: 

1. Solve the equations of motion using the measured input, initial state, and the initial parameter estimate 
6(j to compute the model states and outputs. 

2. Compute the model residual and estimate the noise covariance matrix 


N 


R = diag 


Z=1 


( 5 ) 


3. Compute the output-error cost function, as well as the output sensitivity, cost gradient, and information 
matrices 

S(t) = {^T^^}, for a = 1, 2, ... ,rio and &= 1, 2, ... ,rzp (6) 

f) T 

^ (7) 

i—1 

N 

M = ^S(i)^R-iS(t) (8) 

i=l 


4. 


Update the parameter estimates 
Cramer-Rao bounds 


using the modified Newton- Raphson method^’ ^ 


e = 


00 - M-1 


dJ 


and compute the 


(9) 


cov(0) > M ^ 


(10) 


5. Repeat until the solution converges 

6. Correct parameter uncertainties for colored residuals.^’ 


B. Frequency Domain Estimation 

Output error can be reformulated to use frequency domain data.^’ Measurements are converted into 
the frequency domain by removing the bias and trend, to avoid low-frequency leakage, and then applying a 
high-accuracy chirp-z Fourier transform.®’ This step restricts the equations of motion to linear systems 

x(t) = A(0)x(t) -t B(0)u(t), x(to) = 0 

y{t) = C(0)x(t) -pD(0)u(t) (11) 


3 of 16 


American Institute of Aeronautics and Astronautics 



where the states, control inputs, and model outputs represent perturbations from their respective trim values. 
The system output in the frequency domain is 

y(joj) = [C{jwl - A)-iB + D] u(jw) (12) 

and the cost function to be minimized is 

Nf 

J{e) = Nf^v{jujk)'^S-^v{jujk) (13) 

k=l 

where is the power spectral density of the noise and Nf is the number of frequencies. The output-error 
algorithm then proceeds as normal. 

The most restrictive drawback of using output error in the frequency domain is that only linear systems 
may be considered. Because linear equations describe perturbation dynamics and because of the need to 
detrend the time-domain data prior to Fourier transformation (to avoid leakage), bias parameters cannot be 
estimated in the frequency domain. In addition, practical experience using output error with frequency do- 
main data has exhibited less robustness to model parameter starting values and has required more iterations 
for convergence. 

Offsetting these disadvantages, there are several advantages to using frequency domain data. The Fourier 
transform can be applied at arbitrary frequencies, so that only the data information relevant to the dynam- 
ics of interest is transformed. This typically results in many fewer data points than in the original time 
series, yielding faster computations. Because only the frequency content of interest is transformed, colored 
residual corrections to the uncertainties are unnecessary. In the frequency domain, integration is replaced by 
multiplication, which is computationally faster and makes identifying an unstable dynamic system no more 
difficult than identifying a stable dynamic system. These benefits make the use of frequency domain data 
viable for real-time implementation of output error. 


C. Output Sensitivities 

The output sensitivities reflect how the model parameters change the model outputs, and strongly affect 
the convergence rate and solution. A general method to compute the sensitivities is to use finite differences. 
This method can be applied to arbitrary systems with no change in software. The detriment of doing this is 
that a large number of time-consuming integrations are performed when using time domain data. 

An alternative method is to evaluate analytically-derived sensitivity equations. Although these equations 
are difficult to derive for a general system, for a linear system they are 


d 

dt 


9x\ (9x 5 A 9B 

^ ^ ^ 

89 89 ^ 86^^ 89^ 


9x(to) 

89 


= 0 


The corresponding representation in the frequency domain is 

dy ^ ... 5B A fdC ... dT> .. . 

- = C (M - A) f — x(ju.) + —niju.) ] + f — x(jc^) -f —uijco) 


(14) 


(15) 


This approach involves additional effort to derive and program the sensitivity equations. This effort could 
be considerable when applying to nonlinear systems or when several candidate model structures are being 
evaluated. However, the sensitivities are computed approximately twice as fast using the supplied equations 
rather than numerical finite differences. Additionally, when using noisy data, supplying the model structure 
of the sensitivity equations has a smoothing effect on the sensitivities and therefore helps the nonlinear 
optimization used in output error parameter estimation. These characteristics are well-suited for real-time 
implementation. 


D. Parameter Initialization 

Good starting values of the model parameters are required for timely convergence; otherwise, the algo- 
rithm could use an excessive number of iterations, converge on an incorrect solution, or diverge. Sources for 
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starting values may include previous experiments, wind tunnel tests, an equation-error analysis, or experi- 
ence. 

Denery presented a method in which the traditional output-error algorithm can be initialized at arbi- 
trary values and still obtain quick convergence.^® State histories obtained from measurements, instead of 
integration of the equations of motion, are substituted into the analytically-derived sensitivity equations. 
The parameter estimates then quickly converge to a neighborhood near the solution. 

In addition to adding code complexity, this method requires the sensitivity equations to be analytically 
derived and the state vector to either be measured directly or reconstructed from measurements. Any noise 
or error on the substituted measurements will bias the parameter estimates. However, using this method 
for the first iteration and then switching back to integrating the equations of motion for the remaining 
iterations removes the requirement for good starting values. Practical experience showed that the method 
did not perform as well when using frequency domain data. This is expected to be due to the fact that only 
Fourier transforms are matched, which is the same reason unstable aircraft may be identified using frequency 
domain data with no additional effort. 

E. Sequential Estimation 

Output error can provide real-time estimates by operating sequentially on recent measurements. Time- 
domain data is divided into segments of duration r. While the data system is collecting the next segment 
of data, output error is used to analyze the current segment of data. Parameter estimates and uncertainties 
from the analysis on the previous segment of data is used as prior knowledge for the parameters in a Bayesian 
approach. The output-error cost function, cost gradient, information matrix, and Newton-Raphson updates 
to the parameter estimates are formulated as® 

1 N 1 

J(0) = 2 E + 2 ^p" - (^p) (16) 

F)J ^ 

— = (17) 

i=l 

N 

M = (18) 

0 = 00 - (M + (00 - 0p)) (19) 

where 0p and Sp are the parameter estimates and covariance matrices from the previous analysis. In this 
manner, information from the past is retained. Data forgetting techniques could also be applied here to 
discard past data and more quickly adapt to changing dynamics. 

In the frequency domain, the output data can be transformed using a recursive Fourier transform®’ for 
each frequency w, 

z(jw)j = z(jo;)i_i -f z(i)e“'^“'^‘Af (20) 

and similarly for the input data, where At is the sampling period. An output-error analysis is applied 
periodically using the Fourier transform data in the frequency domain. Parameter estimates from the 
previous analysis are used as starting values, but the previous uncertainty is not used because Fourier 
transforms are cumulative and this would more heavily weight past information. 

Applying output error sequentially yields parameter and uncertainty estimates while flight tests are being 
conducted. There is a trade-off in choosing the size of the data segments. For large segments, sequential 
estimation approaches the traditional post-flight method. For short segments, the information content in each 
data segment is reduced, as well as the amount of time afforded for the output-error analysis, which could 
lead to bad starting parameters and initial conditions, or divergence of the estimator. A good compromise 
is a length that both has enough information content and that provides updates frequently. In this paper, 
the time length r is determined using trial and error with simulation and flight test data. 
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III. Materials 


A. Flight Test Aircraft and Nonlinear Flight Simulation 


The test aircraft used is known as the T-2, which is shown in Figure 1 and described in Table 1. It is a 
5.5% dynamically-scaled version of a generic transport airplane that is used in the NASA Aviation Safety 
Program flight test capability called AirSTAR. Typical flight operations begin with a safety pilot taking the 
aircraft to altitude. Control is then transferred to a research pilot, who flies it from a mobile control room 
using synthetic vision drawn from telemetry data and a database of the local terrain. Control deflection 
commands are generated from the research pilot and the ground-based flight control system. Computer- 
generated excitation inputs for system identification are added to these commands just before the position 
and rate limiters at the actuators. 



Figure 1. T-2 subscale model (credit: NASA Langley Research Center) 


Table 1. GTM mass and geometry values 


Parameter 

Symbol 

value 

Units 

longitudinal center of mass position 

^cm 

0.350 

%c 

mean aerodynamic chord 

C 

0.915 

ft 

wing span 

b 

6.849 

ft 

wing reference area 

S 

5.902 


mass 

m 

1.542 

slug 

roll inertia 

^xx 

1.327 

slug-ft^ 

pitch inertia 

^yy 

4.254 

slug-ft^ 

yaw inertia 

Izz 

5.454 

slug-ft^ 

roll/yaw coupling inertia 

^xz 

0.120 

slug-ft^ 


The aircraft has 16 independent control surfaces. For the present study, the controls were deflected 
in the conventional manner as elevator, aileron, and rudder. Measurements include body-axis translational 
accelerations and angular rates, estimated Euler angles, differential GPS position and velocity, air flow angles 
from vanes mounted on wingtip booms, static and dynamic pressure, fuel flow rate, and control surface 
deflections. First-order anti-aliasing Alters with cut-off frequencies at 16 Hz were used before measurement 
digitization. Data were telemetered to the ground station at 200 Hz, and downsampled to 50 Hz for real-time 
modeling. For additional information about the aircraft and flight operations, see Refs. [19,20]. 

A wind tunnel test article with similar mold lines was studied extensively in the NASA Langley 14- by 22- 
Foot Subsonic tunnel. Static and forced-oscillation data were used to build a high-fidelity nonlinear flight 
simulation called the Generic Transport Model (GTM). A modified version of this simulation, containing 
a smoothed aerodynamic database generated using orthogonal functions, was employed to simulate flight 
data. This simulation was written in MATLAB®and can be run in batch mode with specified inputs or 
flown in real time with a pilot stick and visualization software. 
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B. Flight Maneuver 


The results in this paper use flight test data from a single maneuver that was flown using the T-2 aircraft. 
The measured control surface deflections from this data were supplied as input to the GTM simulation to 
also produce simulated data for that flight. The reference condition is straight and level flight with true 
airspeed 139 ft/s, angle of attack 4.0 deg, and altitude 1226 ft. 

Orthogonal phase-optimized multisines, an input designed and used extensively at NASA Langley Re- 
search Center,®’ were employed to excite the aircraft for system identification. These inputs implement 
simultaneous and uncorrelated excitation on each input over a frequency range of interest and in a manner 
that keeps the aircraft within small perturbations of the trim condition. It is an efficient technique that has 
led to successful modeling results in a variety of unusual flight conditions including hypersonic, high sideslip, 
stall and post-stall flight, as well as normal flight conditions.^®’ 

The measured flight test data from the T-2 aircraft, and also the simulation data from the GTM, are 
shown in Figure 2. Only the longitudinal measurements are shown here, but the lateral and direction axes 
of the aircraft are also being excited in a similar but uncorrelated manner. The elevator input contains a 10 
s excitation with frequencies between 0.2 to 2.2 Hz in 0.3 Hz increments. The input also contains additional 
contributions from the pilot and flight control system. While the responses occur at similar frequencies, 
there are large discrepancies between the flight test and simulation data due to disturbances, differences 
in the test articles, and under-prediction of the damping values in the simulation model. A global Fourier 
smoothing technique® was applied to the flight data to estimate noise covariances, which were then used for 
the simulation data. 




0 2 4 6 8 10 12 


t [s] 


Figure 2. Longitudinal T-2 flight test data and GTM simulation data 


C. Software Development 

SIDPAC is a collection of computer programs written in MATLAB®^® at NASA Langley Research 
Center. This package is bundled with Ref. [5] and contains numerous routines for system identification 
including experiment design, data conditioning, data compatibility analysis, model structure determination, 
parameter estimation, and real-time parameter estimation. Specifically, output-error parameter estimation 
routines using both time and frequency domain data are included in SIDPAC. These codes have evolved 
from previous NASA software,^ new algorithms,®’^®’®® and practical experience applying these methods to 
flight test campaigns at NASA Langley Research Center. 
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The SIDPAC output-error codes were modified to perform the techniques in Section II not already 
included in the codes. Output sensitivities can be computed using numerical forward or central finite- 
differences, or now evaluated using supplied equations. The Denery technique for substituting measured or 
reconstructed state values can be used, and the code can switch back to integrating the state equations after 
a specified number of iterations. The code was already capable of sequential estimation, but timer functions 
were added to limit the execution time during real-time analysis. 

Although real-time operation was achieved, there are several options for increasing code execution speed. 
Results presented herein were performed on a 2.66 GHz laptop. The code could be run on computers 
with more and faster processors. The code could be further streamlined for computational efficiency, and 
parallel computing could be employed to make efficient use of multiple processors. Functions could be coded 
and compiled as MEX-files for increased speed. Alternatively, the entire program could be rewritten in a 
compiled language such as FORTRAN or C for faster speeds. These options were not explored here because 
the MATLAB®version of the code is popular, portable, and can be rapidly changed. Another method for 
quicker convergence times is to change convergence metrics. Currently, there are four metrics that check for 
convergence of residuals, model parameter estimates, cost, and cost sensitivity to within specific threshold 
values. These values are currently set conservatively, and may be relaxed to achieve convergence in many 
fewer iterations but with nearly the same solution. 


IV. Results 


A. Batch Estimation of Short Period Mode 

Output error was first applied to the GTM simulation data in a post-flight analysis to evaluate compu- 
tation time, accuracy, and overall usefulness of the traditional method and its variations. Model parameters 
in the short period approximation 


a 

q 


gS 
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(21) 


were estimated, where flight variables inside the matrices represent the trim values at the start of the 
maneuver. The input vector is augmented with a “1” so that bias terms may be estimated. The model 
parameter vector is 


e = 


Cl^ 


Cl, Cls, 


C„ 


Cr, 


Crr 


07 08 



(22) 


where the last three parameters contain aerodynamic bias terms and trim variables. These bias terms are 
removed from the equation and perturbation states are used when parameter estimation is performed with 
frequency domain data. 

Each analysis was repeated 200 times with unique noise sequences. The initial parameter vector was 
chosen as 6q = [1, 1, 1, — 1, — 1, — 1, 0, 0, 0]^ to start the stability and control parameters with the correct 
sign. Typical fits to the measured outputs are shown in the time domain and the frequency domain in 
Figure 3. In general, the matches were very close, having coefficients of determination above 95% and Theil 
inequality coefficients below 0.06. The slight mismatches between the model and the data may be due to 
the excitation of other modes and small nonlinearities in the data. In all the cases good fits were obtained 
and were in statistical agreement with the numerical linearization solution. 

Parameter estimation results are presented in Table 2, which shows estimated parameters and standard 
errors for and computer execution time, number of iterations until convergence, and time that each 
iteration used. All estimation results were in statistical agreement with the finite difference solutions. Again, 
bias parameters could not be estimated using frequency domain data. Cases 1-4 were performed in the time 
domain, denoted “time,” whereas cases 5-6 were in the frequency domain, denoted “freq.” The output 
sensitivities were computed using finite differences in cases 1 and 5, denoted “numerical,” and the provided 
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(a) time domain analysis (b) frequency domain analysis 

Figure 3. Representative post-flight analysis using short period model 


sensitivity equations were used in the other cases, denoted “analytical.” In cases 2 and 6, the equations of 
motion were integrated to evolve the state vector, denoted “integration.” In case 3, measurements of the 
state vector were substituted for the entire analysis, denoted “substitution,” whereas in case 4 this was done 
for the first iteration only, denoted “switching.” Results for substituting measured state histories in the 
frequency domain data are not shown because the estimates diverged. 

Results for cases 1, 2, 5, and 6 show that computing output sensitivities using finite differences and the 
sensitivity equations yielded the same results and used the same number of iterations. However, using the 
sensitivity equations used about half as much computation time as finite differences. Case 3 used more 
iterations and more computation time than Case 2, and resulted in biased estimates. However, the results 
for case 4 show that substituting state histories for only the first iteration results in the same answer, but in 
fewer iterations and with less computation time. Cases 5 and 6 show that while more iterations were needed 
for convergence in the frequency domain, they were about twice as fast as using time domain data. The 
error bounds here are also lower. 


Table 2. Results of output-error variations using simulation data 


case 

no. 

domain 

sensitivity 

calculation 


07 ± O' 

time [s] 

no. of 
iterations 

iteration 
time [s] 

1 

time 

numerical / integration 

-h4.943± 0.177 

-hO.227 ± 0.008 

2.94 

21 

0.140 

2 

time 

analytical / integration 

-h4.943± 0.177 

-hO.227 ± 0.008 

1.81 

21 

0.086 

3 

time 

analytical / substitution 

-h5.024± 0.192 

-kO.237 ± 0.009 

2.18 

26 

0.084 

4 

time 

analytical / switching 

-h4.943± 0.177 

-hO.227 ± 0.008 

1.66 

19 

0.087 

5 

freq 

numerical / integration 

-k4.958 ± 0.061 

- 

7.60 

90 

0.084 

6 

freq 

analytical / integration 

-k4.958 ± 0.061 

- 

3.27 

90 

0.036 


To further illustrate these results, a typical history of the estimate is shown in Figure 4 for cases 
1-4. Uncertainties are not shown, for clarity. Cases 1-2 follow the same path and converge on the same 
solution. In case 3, the solution quickly approaches the true answer, but overshoots and asymptotes. Case 
4 initially follows the path of case 3, but then switches to evaluating the sensitivity equations, which allows 
it to converge on the same answer as in cases 1-2, but faster and with fewer iterations. 

Based on these results, the following are recommended for real-time implementation of output error: 
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Figure 4. Comparison of estimated short period parameter Clo, evolution 


• Use frequency domain data unless estimating bias terms or the dynamics are nonlinear 

• If practical, provide the sensitivity equations 

• If sensitivity equations are provided, state vector measured or reconstructed, and time domain data is 
used, employ Denery’s method for the hrst iteration to expedite convergence 

B. Real-Time Estimation of Short Period Mode 

In this section, output error is applied to flight data in a real-time manner to estimate the short period 
model. This was done in both the time and frequency domains, using the recommendations of the previous 
section. The GTM simulation was used to compute starting values for the model parameters using numerical 
linearization. The sensitivity equations were provided, and Denery’s method was used for the first iteration 
with time domain data. Corrections were applied to uncertainty estimates when using time domain data. 

The estimation results using time domain data are shown in Figure 5(a) for four different update rates. 
For clarity, only the Cl^ estimate is shown. The hrst plot shows the post-hight solution, which was initiated 
after all the data were recorded and used 2.4 s of computation for convergence. In the second plot, updates 
are requested after collecting 3 s of data. In each analysis, the solver has enough time and information to 
converge on a result. The results evolve and the error bounds decrease as more information is retained, 
ultimately settling on a solution that is in statistical agreement with the batch solution. In the third plot, 
updates are requested every 2 s. The algorithm was interrupted in the majority of these analyses, although 
the estimates were still relatively close and in statistical agreement with the batch solution. In the fourth 
plot, updates were requested every 0.5 s. A 2 s initial delay was afforded because dynamic data is needed 
for parameter estimation. In this case, the output-error algorithm was interrupted during each analysis, and 
outputs were not able to converge. This led to sporadic parameter estimates having large error bounds. The 
algorithm had enough information to converge for the final estimate at the end of the data collection. 

A similar analysis, shown in Figure 5(b), was repeated using frequency domain data. The first plot shows 
the batch solution for the Cl^ estimate. The estimate is in statistical agreement with the time domain batch 
analysis, but has much smaller error bounds and required 6.1 s of computation for convergence. The following 
three plots show results for when updates are requested every 3 s, 2 s, and 0.5 s, respectively. Because the 
frequency domain analysis requires more time for convergence than the time domain analysis, the algorithm 
is interrupted for every analysis. However, the error bounds are relatively small and the parameter estimates 
never stray far from the batch solution, unlike in the time domain. This is because the Fourier transforms 
are cumulative, so that the parameter estimates continue to improve as time goes on. 

In summary, short period aerodynamic model terms can be estimated in real time using output error. 
Accurate results using time domain data are available every 2 s, while results using frequency domain data 
are ready after 3 s of data. If more complicated models with additional parameters are to be estimated, 
these rates could decrease, and of course the update rates depend on the computer being used. 
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0.5s 




(b) frequency domain analysis 


Figure 5. Real-time estimates of short period model parameter Cl^ 
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C. Real-Time Data Compatibility Analysis 


Aircraft sensors contain errors that create kinematically inconsistent sets of data, which then bias aerody- 
namic modeling results. Figure 6 shows flight test measurements, as well as reconstructed values computed 
using the kinematic equations. The mismatch is due to bias and scale factor errors in the measurements. 
Correcting measurements for these errors is called the data compatibility problem, and can be solved using 
output error.® 



Figure 6. Measured and reconstructed longitudinal flight variables 


The sensor measurements are modeled as 


z{i) = (1 + A)y(i) + b + i^{i) 


(23) 


where A is a scale factor and 6 is a bias. For the data compatibility problem, the model parameters are 

iT 


e = 


Av A^ Ac^ by bj^ bf^ A^ A^^ Xip bfp b 0 b^j 


(24) 


although it is common to use only a subset of these. The state equations are found by substituting the 
accelerometer and gyroscope measurements into the rigid body aircraft equations of motion as® 


it 


0 

r - br bq- q 


u 


-g sind + ga^ - ba^ 

V 

= 

br — r 

0 p-bp 
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+ 

g sin (j) cos 9 -\- gay — ba^ 

w 


_ q-bq 

1 

o 

1 

P, 


w 


g cos <j) cos 9 + gttz — baz 


0 


1 sin <j) tan 9 cos 4> tan 9 


1 

1 

1 

9 

= 

0 cos (j) — sin (j) 


q-bq 

_ 


0 sin(/)/cos0 cos (j)/ cos 9 


r — br 


and the model output equations are® 

V = (1 -I- \y)\/u'^ + V^+'w‘^ + by 

P = (1 -I- A,g) arcsin (^v/ \/ u'^ + + bjs 

a = (1 -I- Aq) arctan (w/u) -I- 6 q, 

p = Xcj,)(l> + bp 

9 = {l + \0)9 + bg 

ip = (1 + Xp,)tp bp) 


(25) 


(26) 


Control surface deflections are not included in this analysis because there is no redundant kinematic infor- 
mation for that data. 
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Data compatibility is a difficult estimation problem. The state equations are nonlinear and require 
computationally expensive and iteration solvers. Good starting values are needed and the algorithm can 
converge on incorrect solutions due to low sensitivities. Sufficient excitation is required to accurately estimate 
the model parameters, e.g. the airspeed scale factor would be difficult to estimate using data from a short- 
period maneuver, where airspeed is practically constant. 

One approach to solving the data compatibility problem in real time is to linearize the system about 
the trim condition. To illustrate, applying partial differentiation to Eqs. 25 and 26, evaluating at the trim 
conditions for straight and level flight, and specializing to the longitudinal dynamics yields the linear system 

Ail 0 —bq —gcosOo Au 1 0 —wq Aux 

Aw = bq 0 —g cos 4>o sin 6q Aw + 01 uq Aa^ 

A9 0 0 0 A9 0 0 cos (j)Q Aq 

’ ad 1 r mo(1 + Av)/Do wo(l + Ay)/Do 0 1 T Aw ‘ 

Aa = -wq{1 + Xo)/V^ uo(l + Aa)/Do^ 0 Aw (27) 

A9 \ [ 0 0 1 + J [ Ad 

Although the scale factors and gyroscope biases can be estimated using this system, the other biases do not 
appear in the equations and cannot be estimated. If data corrected using these equations was then used for 
aerodynamic parameter estimation, only the stability and control derivative estimates would be free of data 
compatibility related errors. 

Another approach is to use the full nonlinear data compatibility equations in Eqs. 25 and 26. A time- 
domain analysis is required because the equations are nonlinear, but every unknown parameter can be 
estimated. In this case, the full data set can be corrected for data compatibility errors for further analysis. 
This approach is presented next, using flight data from the T-2 aircraft. The output sensitivity equations, 
derived by hand and checked using a computer algebra software called Maxima,^® were supplied to expedite 
the estimation. Denery’s method for state substitution was used for the first iteration to improve the 
convergence. Only the longitudinal equations were used in this case. Biases for the accelerometer and 
gyroscope, and scale factors for the angle of attack and pitch angle were estimated, which is typical. 

Figure 7 shows the model parameters, estimated every 4 s. The parameters and uncertainties eventually 
settle with each segment of data. The final model parameter estimates are small because of high quality 
measurements and calm air during the flight test. The sensor measurements were corrected for each segment 
and are shown in Figure 8. Except for the longitudinal acceleration signal, these traces are almost indistin- 
guishable. Estimation of the short period model parameters using this corrected data resulted in less than 
a 1% change in stability and control parameters and no change in the estimated uncertainties. 

The 4 s update period was found to be a good trade off between accurate results and the update rate 
of the model parameters. The largest computational burden was the integration of the nonlinear equations 
and the uncertainty correction for the colored residuals. One option would be to check the whiteness of 
the residuals and skip this correction, because the uncertainties are only of secondary importance here, and 
because the kinematic equations are well known and can only be corrupted by process noise and unmodeled 
sensor errors. In this case, accurate model parameters could be estimated up to every 0.5 s. 

V. Conclusions 

Output error is a standard method for estimating unknown parameters in dynamic models of aircraft 
in post-flight analysis. This is due to a sound theoretical underpinning, favorable asymptotic properties, 
realistic assumptions consistent with the real flight test situation, and a straightforward approach that has 
been tested and improved upon using numerous flight test cases. However, there are drawbacks that have 
prevented its application to real-time parameter estimation during flight. 

This paper addresses some of those difficulties by exploring variations of the traditional output-error 
approach that make it more suitable for real-time applications. Real-time estimation is becoming increasingly 
important because it helps to reduce time and costs associated with identifying dynamic models from flight 
test data. For instance, immediate modeling feedback can shorten maneuver time durations or help diagnose 
problems while the aircraft is still on condition, instead of after the flight test. 

Experiments using simulation data to estimate short period model parameters led to several suggestions 
for real-time implementation. Frequency domain data should be used unless the dynamics of interest are 
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nonlinear and cannot be linearized, or if bias terms are of interest. If practical, the output sensitivity 
equations should be derived and provided to the estimator to increase computation speed. In this case, 
Denery’s method can be used with time-domain data to expedite convergence from arbitrary starting values 
for the model parameters. Application to flight test data using MATLAB®code and a laptop computer 
demonstrated that short period model estimates could be estimated every 2 s in the time domain, or more 
rapidly in the frequency domain after 3 s of data. A list of recommendations for increasing computational 
speed was provided. Model size (number of states, outputs, and parameters) as well as the capability of the 
computing equipment affect the achievable update rate. 

Application to flight test data also showed that the data compatibility problem could be solved every 4 s, 
or as fast as every 0.5 s if the error bound corrections for colored residuals are skipped. Data compatibility 
is an important and standard step in system identification that ensures a kinematically consistent set of 
data, prior to aerodynamic modeling. In the case presented, the corrections were small and only led to a 1 % 
change in the estimated short period estimates. However, in other cases this is not the case and corrections 
could be much larger. A 4 s update rate may be slower than desired because the aerodynamic estimation can 
update faster, but again the uncertainty correction can be skipped and the code run time can be shortened. 
In addition to correcting data for aerodynamic modeling, solving the data compatibility problem in real time 
can be used for flight operations data checks or for fault detection. If estimated bias or scale factors are 
abnormally large, the flight test engineers can be alerted to a potential fault in the associated measurements. 

There still remain shortcomings of applying output error in real time. A model structure is needed before 
modeling; however in many instances, such as perturbation maneuvers, low-order equivalent system modeling, 
and data compatibility analysis, an adequate model structure can be assumed. Due to assumptions, output 
error is still not equipped to deal with the effects of turbulence and other state disturbances. Bias parameters 
can only be estimated using time domain data because of detrending prior to transforming the data into 
the frequency domain. Output error still is fundamentally an iterative technique that requires more time 
than non-iterative methods. Finally, a shortcoming of Denery’s method using arbitrary starting values of 
the model parameter estimates was identified when using frequency domain data. 

This method provides another approach to real-time dynamic modeling as the aircraft is flying. Bias 
parameters, such as the profile drag, can be estimated directly in real time. The data compatibility analysis 
will provide flight test engineers with real-time updates to sensor measurements, and the ability to check 
for proper sensor function in real time. The data compatibility analysis could be integrated into future 
algorithms for fault detection. This method could also be used for other modeling applications, such as 
fitting transfer functions to real-time frequency response data. Data forgetting techniques could also be 
incorporated to make the algorithm more rapidly adapt to changing dynamics. 
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